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!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! 

I’d like to remind you that TIME 2012 is approaching. Please consider participa- 
tion and don’t forget to submit a lecture and/or workshop. Marina Lepp and Eno 
Tonisson are working on the social program. Fortunately we could win a promi- 
nent lady from Estonia for giving an additional keynote lecture. 

Bitte iiberlegen Sie die Teilnahme an der TIME 2012 und vergessen Sie nicht, einen Vor- 
trag oder Workshop anzumelden. Die erste Einreichung ist bereits eingelangt (aus Austra- 
lien! !). Marina Lepp und Eno Tonisson arbeiten bereits am Begleitprogramm. Erfreuli- 
cherweise konnte auch eine prominente Dame aus Estland fur einen zusatzlichen gewon- 
nen werden. 


Attend TIME2012 in Tartu, Estonia 
July 11-14, 2012 
http://time201 2.ut.ee/ 

The next announcement is for now for our German speaking members: 

Vorerst ergeht diese Einladung an unsere deutschsprachigen Mitglieder: 

Mitgliedschaft bei der ACDCA 

Das ACDCA (Austrian Center for Didactics of Computer Algebra) wurde Ende der 80er 
Jahre gegriindet und ist seither eine der treibenden Krafte fur den Einsatz von CAS in der 
M-Ausbildung. Dieser Verein wurde kurzlich umstrukturiert und bietet alien Interessier- 
ten die Mitgliedschaft und die Gelegenheit zur Mitwirkung an. Bitte informieren Sie sich 
auf der Homepage des ACDCA iiber die Aufgaben und Ziele des ACDCA. 

Info und Beitrittserklarung linden Sie auf 

http://www.acdca.ac.at 

Our DUG member Klaus Korner is asking for a FORTRAN Compiler. 

It needs not to be a fancy one. There are several Freeware Compilers 
offered on the web. Can anybody provide a special recommendation? 
Please write and I will forward your mail to Klaus. 

Many thanks in advance 

Klaus has an interesting proposal for a CAS application. There is not enough 
space in this DNL to write more. I will present his idea in the next issue, 

Josef 
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Liebe DUG-Mitglieder, 

Ich bin sehr froh, dass es mir dieses mal 
gelungen ist, den DNL rechtzeitig fer+ig 
zu stellen. 

Dieser DNL besteht zum gro'Bten Teil aus 
Reaktionen unserer Mitglieder auf fruhe- 
re Beitrage und auf Anfragen zu DERIVE 
und zum TI-Nspire. Fur den Herausgeber 
ist es naturlich sehr erfreulich, dass die 
DNL-Artikel zum Weiterforschen anre- 
gen. 

Stefan Welke lieferte eine tiefgrundige 
Behandlung eines Aufsatzes von Carl Lein- 
bach und mir und unser Freund David 
Halprin aus „Downunder“ erinnert sich an 
seine Beschaftigung mit Flugzeugprofilen. 

Eine Anfrage von Robert Setif traf sich 
mit einer Frage von Rob Gough zum popu- 
laren Thema „Mandelbrotmenge“. Dies war 
auch fur mich Anlass genug, mich mit DE- 
RIVE daran zu versuchen. 

Ich darf berichten, dass die DUG auch 
heuer wieder neue Mitglieder gewinnen 
konnte. Erik von Lantschoot - ein neues 
Mitglied - hat sich auch gleich mit einem 
trickreichen TI-Nspire-Programm zur 
Dreiecksauflosung eingestellt. Vielen Dank 
dafur. 

Piotr Trebisz zeigt die dritte Folge seiner 
Schneckenhtiuser und schlieBlich wiinscht 
Roland Schroder mit einem vierbltittrigen 
Klee den DUG-Mitgliedern ein erfolgrei- 
ches neues Jahr 2012. Diesen Wunschen 
schlieBe ich mich gerne an. 

Schone Beitrage aus Kolumbien (Nelson 
Urrego) und USA (Phil Todd) sind wieder 
eingelangt. 

Beachten Sie bitte die Informationen auf 
der gegenuberliegenden Seite. 

Ich verbleibe mit den besten GruBen und 
Wunschen und mit der Hoffnung auf ein 
Treffen in Tartu (TIME 2012). 


Dear DUG Members, 

I am very happy that I could finish 
DNL#84 in time. 

The major part of this DNL consists of 
reactions of our members on earlier con- 
tributions and on requests to DERIVE and 
TI-NspireCAS. It is very enjoyable for 
the editor that DNL articles inspire for 
further investigations and for providing 
comments. 

Stefan Welke sent a profound treatment 
of an earlier paper presented by Carl lein- 
bach and me. Our friend David Halprin 
from Downunder remembers his work with 
airplane prof iles “some" years ago. 

A request from Robert Setif met a ques- 
tion posed by Rob Gough both dealing with 
the popular issue “Mandelbrot Set". This 
was a welcome occasion for me to make a 
try with DERIVE. 

I can report that the DUG could win again 
new members in 2011. Erik van Lantschoot 
- a new member - presents a tricky TT- 
Nspire-program for solving triangles. 
Many thanks for this. 

Piotr Trebisz shows the 3 rd part of his 
snail shells and finally Roland Schroder 
sends a four-leaf clover wishing all DUG 
members a successful New Year 2012. I 
would like to join him wishing you and your 
families the very best. 

Great contributions from Colombia (Nel- 
son Urrego) and USA (Phil Todd) have 
arrived. 

Please notice the information on the oppo- 
site page. 

I remain with best regards and wishes and 
hope to meet you in Tartu (TIME 2012). 



Download all ZWL-Derive- and Tl-files from 

http : / / www . austromath . at/ dug/ 
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Contributions: 

Please send all contributions to the Editor. 
Non-English speakers are encouraged to 
write their contributions in English to rein- 
force the international touch of the DNL. It 
must be said, though, that non-English 
articles will be warmly welcomed nonethe- 
less. Your contributions will be edited but 
not assessed. By submitting articles the 
author gives his consent for reprinting it in 
the DNL. The more contributions you will 
send, the more lively and richer in contents 
the DERIVE & CAS -TI Newsletter will be. 


Next issue: March 2012 


Preview: Contributions waiting to be published 

Some simulations of Random Experiments, J. Bohm, AUT, Lorenz Kopp, GER 
Wonderful World of Pedal Curves, J. Bohm, AUT 
Tools for 3D-Problems, P. Luke-Rosendahl, GER 
Hill-Encription, J. Bohm, AUT 

Simulating a Graphing Calculator in DERIVE , J. Bohm, AUT 

Do you know this? Cabri & CAS on PC and Handheld, W. Wegscheider, AUT 

An Interesting Problem with a Triangle, Steiner Point, P. Luke-Rosendahl, GER 

Overcoming Branch & Bound by Simulation, J. Bohm, AUT 

Graphics World, Currency Change, P. Charland, CAN 

Cubics, Quartics - Interesting features, T. Koller & J. Bohm, AUT 

Logos of Companies as an Inspiration for Math Teaching 

Exciting Surfaces in the FAZ / Pierre Charland's Graphics Gallery 

BooleanPlots.mth, P. Schofield, UK 

Old traditional examples for a CAS - what's new? J. Bohm, AUT 

Truth Tables on the TI, M. R. Phillips, USA 

Where oh Where is It? (GPS with CAS), C. & P. Leinbach, USA 

Embroidery Patterns, H. Ludwig, GER 

Mandelbrot and Newton with DERIVE , Roman Hasek, CZK 

Tutorials for the NSpireCAS, G. Herweyers, BEL 

Some Projects with Students, R. Schroder, GER 

Dirac Algebra, Clifford Algebra, D. R. Lunsford, USA 

Treating Differential Equations (M. Beaudin, G. Piccard, Ch. Trottier), CAN 

A New Approach to Taylor Series, D. Oertel, GER 

Statistics with Tl-Nspire, G. Herweyers, BEL 

Cesar Multiplication, G. Schodl, AUT 

Henon & Co; Find your very own Strange Attractor, J. Bohm, AUT 
Rational Hooks, J. Lechner, AUT 
Cubus Simus, H. Ludwig, GER 

Simulation of Dynamic Systems with various Tools, J. Bohm, AUT 

Using DERIVE to simulate the basic steps of the Q. Shor’s algorithm, N. Urrego, COL 

Approximating Circular Arcs with Cubic Splines, Ph. Todd, USA 

and others 

Impressum: 

Medieninhaber: DERIVE User Group, A-3042 Wtirmla, D'Lust 1, AUSTRIA 
Richtung: Fachzeitschrift 
Herausgeber: Mag. Josef Bohm 


The DERIVE-NEWSLETTER is the Bulle- 
tin of the DERIVE & CAS- 77 User Group. 
It is published at least four times a year 
with a content of 40 pages minimum. The 
goals of the DNL are to enable the ex- 
change of experiences made with DERIVE, 
TI- CAS and other CAS as well to create a 
group to discuss the possibilities of new 
methodical and didactical manners in 
teaching mathematics. 

Editor: Mag. Josef Bohm 
D'Lust 1, A-3042 Wiinnla 
Austria 

Phone: ++43-(0)660 3136365 

e-mail: nojo.boehm@pgv.at 
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Here p ist he vector of coefficients p= [p, q, r] from (1.3), n= [nl , n2 ] gives the first and 
last exponent in {M nl v 0 ,M" l+l v 0 ,...,M n2 v 0 } , and v= [ z 1 , z 2 , z 3 ] the initial vector. 

#2: g( [0.5, 0.5,0] , [0,7] , [0 , #i, 1 ] ) 

#3: 

[[0,0], [0,1], [1,0], [0,0.5], [0.5, 0.5], [0.5,0.25], [0.25,0.5], [0.5,0.37 
5] ] 


The objective of Carl’s and Josefs article [1] 
was the proof of convergence of (1.1) respec- 
tively. (1.2). We will conclude this introduction 
with two simple propositions about limits of 
recurrences. 


Proposition 1.1 

If the limit z X: := lim z n + 0 of (1.2) exists, then 

n — >oo 

p + q + r = 1 . 

Proof: If this limit exists, then obviously z x := limz n = limz n+1 = lim z n+2 = lim z n+3 + 0 . Thus 

n—> oo n—> oo n — >oo n—> oo 

z x = p-z x +q-z x +r-z x = (p + q + r)z x . Division by z x + 0 gives the result. □ 

As we shall see, p + q + r = 1 is necessary but not sufficient for convergence. The next 
proposition is about boundedness of the recurrence. 

Proposition 1.2 

1.2.1. If under the conditions (1.2) \p\ + \q\ + |r < 1, then the sequence { z „}„ eH Z5 
bounded by M := max{|zj | , |z 2 1 , |z 3 1 } and has at least one accumulation point. 

1.2.2. If under the conditions (1.2) Ipl + U + |/j = Q < 1 , then limz n = 0 . 

II II II x— »oo 

Proof of 1.2.1.: From the triangle inequality follows: 

k + 3 1 ^ \P z n + 2 1 + \ c l z n + 1 1 + \ rz n | < (|^| + 1 ^| + 1 ^|) • max {|z„ +2 1 , | z„ +1 1 , \z „ |} < max {|z„ +2 1 , | z n+1 \,\z„\}=: M n 

An immediate consequence is M n+1 = max{|z n+3 |,|z B+2 |,|z n+1 |} < max{|z n+2 |,|z B+1 |,|z B |} =M n . We 
conclude by induction: | z n | < M n <M X =M. This proves the first part. 

The boundedness implies, that the sequence stays within the closed polycylinder 
D := {(z (1) ,z <2) ,z (3) ) e C 3 | |z (y) | <M,j e {1,2,3}} , which is a compact subset of C 3 . This 

proves the second part. 
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Proof of 1.2.2.: The same reasoning as above gives M n+3 < Q-M n , so by induction 
\z 3n+k I < M 3n+k < Q " -M x =Q n -M for k e {0,1,2} . Since 0 < Q < 1 , 1.2.2. is proved. □ 


2 

The third proposition is about the convergence of the sequence of quotients < > . As we 

1 n J 

know from the Fibonacci numbers, the sequence of quotients may converge, even if the se- 
quence itself does not. 

Proposition 1.3 

If the sequence of quotients of (1.2) converges, and if the limit q x := lim ^ is not 

n— >oo n 

equal to zero, then q y satisfies the algebraic equation q j - r ■ qf -q -q x — p = 0 . 


Proof: Assume, the limit exists, then clearly q r := lim — = lim = lim ^ follows, because 

n — >oo n n — >oo W+1 n—> oo M+2 

almost all numbers of the recurrence are not equal to zero. We now divide (1.3) by z n and 


obtain the equivalence = p + q + r ■ - !} ^ L 


p + q ^±L + r ^.^±I 


Passing to the limit, the equation becomes qj = r ■ qj + q ■ q y + p . This is equivalent to the 
statement of the proposition. □ 


Remark: If we choose z 0 = 0, z x = 1, z 2 = 1 , this three-step recurrence produces with p = 1 , q = 2 , 
and r = 0 , as well as with p = 0 , q = 1 , and r = 1 and with any linear combination of these two 
parameter vectors (p, q, r) = a- (1, 2, 0) + (1 - a) ■ (0, 1, 1) with a e R the Fibonacci numbers. We 
give an example: 


#4: 

W ~ 0.7* 

[i. 

2, 

0] 

+ 0.3- 

[ 0 , 

1, 

1] 

#5: 

g(v, [0, 

10] 

, [ 

0 , 

1, 1] ) 

1 



#6: 

[ 0 , 

1, 

1, 

2, 

3, 5, 

8, 

13, 

21, 34, 55] 


(2) The Binet - Representation 

We compute the determinant and the characteristic polynomial of the matrix M from (1.5): 

#7: M == [0, 1, 0; 0, 0, l;p, q, r] 

#8: DET(M) 

#9 : p 

#10: CHARPOLY (M, A) 

3 2 

-A+r-A +q-A+p 


# 11 : 






Stefan Welke: Using Linear Algebra to Explain . . . 


D-N-L#84 



We notice, that the possible limit of the sequence of quotients from proposition 1.3 must be a root of 
the characteristic polynomial of M . 

Now let M have three distinct roots, X^,X 2 , A , possibly one equal to zero. These roots are the eigen- 
values of M . Let X be any eigenvalue, then clearly U = p + qX + rX 2 . This explains the following 
identities: 

f° 1 0 YM f 2 If 1 ! f 1 ! 

(0.6) 001 X = X 2 = X 2 =X- X 

y p q r jyX 2 ) [^p + q ■ X + r ■ X 2 J yX* ) 

Therefore the corresponding eigenvectors are, up to a factor: 

f 1 ! f 1 ! f 1 ! 

(0.7) e x = X I ,e 2 = X 2 ,e 3 = X 3 

A j A j A j 

We assume, that all eigenvalues are distinct, then Vandermonde’s determinant gives 

f 1 1 1 ! 

(0.8) det X 1 X 2 X 3 =(X l -X 2 )(X 2 -X 3 )(X 3 -X 1 )*0 . 

A'AA 2 , 

This result proves the linear independence of the three eigenvectors (0.7) for distinct roots, 
X l ,X 2 ,X 3 , and we can represent every initialvector v 0 e C 3 as a linear combination of (0.7), 

because the vectors (0.7) form a basis of C 3 : 

(0.9) v 0 = a x e x + a 2 e 2 + a 2 e 2 , a k e C 3 

We call (0.9) the ih'net-decomposition of v 0 . If we now calculate (1.5) with this decomposition, then 
we obtain by linearity: Mv 0 = M(a l e l + a 2 e 2 + a 2 e 2 ) = a 1 X l e l + a 2 X 2 e 1 + a 2 X 3 e 2 and hence 

(0.10) v n =M"v 0 =M n (a l e l +a 2 e 2 +a 2 e 2 ) = a l X l "e l +a 2 X 2 n e 2 +a 2 X 3 e 2 

At this point, after passing to the first coordinates of the vectors v 0 and v n , the proof of the 
first part of the following proposition is completed: 

Proposition 2.1 

Assume, that the characteristic polynomial c(X) = X 3 - rX 1 -qX- p has three distinct 
roots X l ,X 2 ,X 3 , and a x ,a 2 ,a 3 are the coefficients of the Binet-decomposition of the 
initialvector, then the linear recurrence (1.3) has the following representation: 

(0. 1 1) z n = a x Xf + a 2 X 2 " + a 2 X 3 , neN 


This representation is unique. 
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The uniqueness of this representation is a consequence of the uniqueness of the eigenvalues of 
a matrix and of the uniqueness of the representation of a vector with respect to a given basis 
of the vector space. 

We call (2.6) the Binet- formula, or the 5me/-representation of the recurrence. Let us look at 
our example with p = q = \, r = 0 , and the initial values z 0 =0,z 1 = /, z 2 = 1 . The roots of the 
characteristic equation and the coefficients of the /ime/-decomposition are: 

#12: InputMode : = Word 

#13: SOLVE (A A 3 - 0.5-A - 0.5 = 0, A) 

#14: A = -0.5 - 0.5-i OR A = -0.5 + 0.5-i OR A = 1 

#15: SOLVE ( [ 0 , 1, i] = al-[l,-0.5 - 0.5-i, (-0.5 - 0.5-i) A 2] + 

+ a2-[l,-0.5 + 0.5-i, (-0.5 + 0.5-i) A 2] + a3-[l, 1, 1], [al, a2, a3] ) 

#16: al = 0.4 + 0.2-i AND a2 = -0.8 - 0.6-i AND a3 = 0.4 + 0.4 -i 

Thus the Binet - representation of Carl’s and Josefs recurrence is 

(0. 12) z n := (0.4 + 0.2/)(-0.5 - 0.5/)" + (-0.8 - 0.6/)((-0.5 + 0.5/)" + 0.4 + 0.4/ . 

The limit of this sequence is easily found, because |— 0.5 - 0.5/| = |— 0.5 + 0.5/| = ^ < 1 . We get: 
lim z n := (0.4 + 0.2/) • lim(-0.5 - 0.5/)" + (-0.8 - 0.6/) • lim(-0.5 + 0.5/)" + 0.4 + 0.4/ = 0.4 + 0.4/ 

n — >oo n—>cc n—> oo 

A close examination of the 5me/-representation (0.1 1) tells us that 

• the convergence or divergence of the recurrence depends on the eigenvalues of the 
characteristic polynomial, hence on the parameters p, q, r in the definition of the re- 
currence, and on the initial values of the sequence. 

• there exits an interpolating function b(t ) := a x / 1/ + a 2 A 2 + a 2 /J ,t e I , which is the su- 
perposition of at most three spirals, if |/L y + 1 , and/or circles, if \A k | = 1 . 

We will give an illustrative example: 

Choose three eigenvalues ij =-l,A 2 = -1 + i,A 3 = -0.5 - 0.5/ , then DERIVE gives: 

#17: EXPAND ((z + l)-(z + 1 - i)-(z + 0.5 + 0.5-i), z) 

#18: z A 3 + 2.5-z A 2 + 2.5-z + 1 - i-(0.5-z A 2 + 0.5-z) 

So the parameters are p = -l,q = r = -2.5 + 0.5/ . The next function spiral l (w, t, v) com- 
putes the parametric equation of the spiral with the parameter vector w= [p, q, r] , the variable 
t , and the initialvector v= [ z o , z l , z 2 ] : 

#19: spirall (w, t, v) — 

Prog 


e == EIGENVALUES ( [ 

[0, 

l. 

0] , 

[0, 0, 1], 

w] ) 

w == VECTOR ( SUBST ( 

[l; 

s; 

s A 2] 

, s, eik). 

k, 1, 3 

z == SOLUTIONS ( [a. 

P, 

Y] 

• w = 

[ v | 1 ; v|2; 

vj3] , 


[a, (3, y] ) -VECTOR (If (k + 0, k A t, 0, 0), k, e) 
[RE ( z ) , IM ( z ) ] ' il 


#20: spirall ( [-1, -2.5 + 0.5-i, -2.5 + 0.5-i], t, [0, 1, i] ) 
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The ihheZ-decomposition is: 


f (U 


r r 


f 1 N 


f 1 ^ 

1 

JJ 

= 3 

-i 

“(s + sO 

-1 + / 
l ~ 2i > 

1 

“>|5 

1 

1 

2 2 
i 

K 2 J 


The eigenvalue X 2 with |-l + z| = >/2 > 1 is 
responsible forlimh(t) = oo , i.e. that the spiral 

t — >oo 

tends to infinity. 


We change now the set of initial values, [1,-0. 5-0. 5 i , 0 . 5 i ] in the first case: 


#22: spirall ( [-1, -2.5 + 0.5-i, -2.5 + 0.5-i], t, [1, -0.5 - 0.5-i, 0 . 5 - i ] ) 
#23 : [2 A ( -t/2 ) -COS (3-n-t/4) , -2 A ( -t/2 ) -SIN (3-n-t/4) ] 


and [ l , - 1 , l ] in the second case: 


#24: spirall ( [-1, -2.5 + 0.5-i, -2.5 + 0.5-i], t, [1, -1, 1]) 



#25 : [COS (n-t) , SIN (n-t) ] 

In the first case the ZhneZ-dccomposition is 

f 1 

v o =1 ‘ -J-J 

V 2 j 


so b(t) := , is the Binet- 

representation, which is the equation of a 
clockwise winding spiral with limit point 
0 . In the second case we find 


r n 


=i- 


and b(t) = (-1)' = e ,7It , 


which is the parametric form of the unit 
circle in C . 


A last example in this instance is a combination of the last two cases with 


r n 


r i 

'N 


( 6 h 

5 

-i 

-i- 

i 

i 


— 1 — -L / 
10 2 1 

2 

2 


vU 


i 

\ 2 

y 


! + !/ 

5 ^ 2 l J 
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The corresponding curve winds 
clockwise around a circle with ra- 
dius } . If we had added any multi- 
ple of the third eigenvector, how 
ever small in absolute value, we 
would see instead something like a 
spiral tending to infinity. 

These examples demonstrate the 
impact of the initial values. 

The following DERIVE - function is 
part of the function spiral l and 
computes the coefficients of the 
Z?//ze/-dccomposition of a initial 
vector. 


#27: binet_decomp (p, v) — 

Prog 

e ~ EIGENVALUES ([ [0, 1, 0], [0, 0, 1], p] ) 

w ~ VECTOR (SUBST ( [1; s; s A 2], s, e;k), k, 1, 3) 

[e, SOLUTIONS ( [a, (3, y]-w = [vj.1; vj.2; vj.3], [a, (3, y] ) ] 


#28: binet_decomp ( [-1, -2.5 + 0.5-i, -2.5 + 0.5-i], [0, i, 1]) 

#29: [[-1, -1 + i, -0.5 - 0.5-i], [[3, -0.6 - 0.8-i, -2.4 + 0.8-i]]] 


(3) Constructing Recurrences to Given Figures 

We are in the position to explain the occurrence of the ellipses in Josefs example [l,p. 24], 
because we have a theory, that can predict the resulting phenomena. Consider the case of 
three distinct eigenvalues \\ | < 1, \A 2 1 = \A 3 1 = 1, A 2 ^ I and A 2 =A i , where the coefficients of the 

5me/-decomposition satisfy the conditions a 2 -a 3 ^0 and \a 2 \ ^ \a 3 \ . Then the curve given by 

c(t ) := a 2 ■ A 2 + a 3 ■ A[ is an ellipse. The reason is simply, that C is a two-dimensional vector 
space over R. , and that a real linear map / : C — > C has the form / (z) = a ■ z + J3 • z with 
a,j3 e<C . The image of a circle is in general an ellipse, and it is evident, that we have 
c(t ) = /(A' ) , with f(z) := a 2 ■ z + a 3 ■ z . The conditions on a 2 , a 3 guarantee, that we obtain 
an ellipse, which is not a circle. The following example illustrates this case: 

Let \ = 0.5, A 2 = -0.8 + 0.6/ and A 3 = -0.8 - 0.6/ be the roots of the characteristic polyno- 
mial. Then the parameters of the recurrence are p = 0.5,7/ = -0.2,/- = -1.1 . We choose three 
distinct initial vectors and calculate the corresponding //meZ-dccompositions: 



20 + 32 i 

63 - 34z 

57 + 226; 

II 

II 

f 1 ) 
-0.8 - 0.2; 


C l 

C 2 

c 3? v 2 


41 

246 

246 

v 0.28 + 0.32;'y 
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Now two cases can happen: 


with 



f 1 'l 


f 1 1 


= 

0.5 


-0.8 + 0.6/ 

> e 3 = 




^0.28-0.96;, 



The corresponding curves are red 
for v, , blue for v 2 , and black for 


The recurrence (1.3) produces a 
sequence of points {z 0 ,z 1 ,z 2 ,...} , 

and we can draw the correspond- 
ing sequence of line segments, 
each segment connecting two 
subsequent points z n and z n+1 . If 
the curve given by b(t) is a circle, 
then all segments have equal 
length. 


(1) The sequence is periodic, then the sequence of line segments looks like a regular polygon. 
This happens, if the corresponding eigenvalue is of the form X = e a n ' ‘ , a e Q \ {0} . 

(2) The sequence is not periodic and the points are a dense subset of the circle. In this case the 
corresponding eigenvalue is of the form X = e a ' n i , a e M \ Q . Then the envelope of all line 
segments is a smaller circle with the same centre as the generating circle. 


In the case of two conjugate eigenvalues of absolute value 1 we thus have as envelopes of the 
line segments either a polygon with vertices on an ellipse (1), or an ellipse (2). We will inves- 
tigate Josefs example from [l,p. 24]. The first observation is, that Rudiger Baumann’s ap- 
proach is almost equal to our first function g (p, n, v) , but he uses the iterates function in- 
stead of matrix powers. We recognize, that in our notation p = 0.95, c/ = -0. 1 and r = 0 in the 
first example of [ 1 ,p. 24] and p = -0.95, q = -0. 1 and r = 0 in the second. Now we look at 
the roots of the characteristic polynomial: 


#30: APPROX (SOLUTIONS (X A 3 + 0.1 -X + 0.95 = 0, X)) 

#31: [-0. 94915332 40, 0.4745766 62 0+0. 8807207416 ■ i, 0.4745766 62 0-0. 8807207416 • i] 

#32: APPROX (ABS (#32 SUB 2)) 

#33: 1.000445916 

#34: APPROX (SOLUTIONS (X A 3 + 0.1 -X - 0.95 = 0, X)) 

#35 : [0. 94 915332 40,-0.47 457 6 6620 + 0 . 8807207416 • i, -0.4745766 62 0-0. 8807207416 • i] 

The real eigenvalue has absolute value smaller than 1, so its contribution to the Binet- 
representation is negligible for large values of n . In both cases the other two eigenvalues are 
conjugate to each other and of absolute value nearly 1. 
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We finish this section with two examples. The first is a periodic recurrence. Set /t, = 0 , 
a 2 = e l ' n/1 , and A 3 = e '' n/1 . We compute the parameters p,q,r with fl and then draw the 
picture with polygonl: 


#36: fl(x, y, z) — [x-y-z, - (x-y + y-z + z-x), x + y + z] 

#37 : fl (0,EXP(pi-i /7 ) , EXP (-pi • i/7 ) 

#38 : [0, -1, 1.801937735] 

#39: pointseql (p, n, v) — VECTOR (spirall (p, k, v) , k, n SUB 1 , n SUB2 ) 



The next example demonstrates that we can 
even achieve cycloids and periodic polygons 
with vertices on a cycloid. We use the Binet- 
representation to construct a recurrence with 
the desired property. 

We want to construct a cycloid with three 
cusps on a circle. A possible parametric repre- 
sentation in complex notation is: 
b(t) = 4 • e iat + e Mat ,ael\{0} 

To construct a recurrence with the above Binet- 
representation, choose eigenvalues 
4=1,^= e'“= 0.8 + 0.6/, A 3 = e 4, a = A , 4 . We 
perform the following computations: 


#40: 1 : =0 . 8+0 . 6 • i 
#41: fl (1, a, a A 4) 

#42: f2 (x, n, m) == n-[l, x, x A 2 ] + [1, x A m , x A (2m) ] 

#43: f2 (a, 4, 4) 

#44: [0. 151243161 5+0. 9884965 887 -i, -1. 5332 15 64 1-2. 175105228-i, 2. 381 97248+ 



#45: spirall (#42, t, #44) 

Note, that the function f 2 (x, n, m) computes the 
initial vector v 0 = 4 ■ < 
the figure to the left. 


initial vector v 0 =4-e 2 + e 3 . The result looks like 


If we want to construct a recurrence with a peri- 
odic polygon with vertices on a cycloid, we need 
to choose eigenvalues, which are roots of unity, 
e.g.: 

\ — 1, a 2 — £ ? /l^ — 6 — Ai 

We omit the calculations and just present the result 
in the next figure: 
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The previous examples suggest the question, 
which figures are constructible with recur- 
rences. So far, we have been able to construct 
spirals, circles, ellipses, and cycloids as inter- 
polations of ///net-representations. The result- 
ing sequences are periodic or infinite. 

We proceed to construct recurrences, where the 
corresponding 5/«e/-reprcscntation leads to a 
parametric curve, which is either a parabola or 
a hyperbola. 

Let us begin with a parabola. To do so, we 
choose three real eigenvalues A, > 0, 2. = a; , 
and T, = 1 . Now choose the initial vector 




f 1 > 


f n 

(0.13) 

II 

_ r*-) 

+ 

<N 

II 

O 

A 

+ i 

A 2 





l^ 4 J 


Then the interpolating function is b(l ) = A + iA 2 = A + / ( A j . The resulting parametric curve is 
part of a parabola as the next figure shows. 

Take A = 2, A = 4 . Then (z- l)(z -2)(z-4) = z 3 -7 z 2 +14z - 8, so p = $>,q = -14 , r = 7 . We 



rn 


'n 


' 1 + / ^ 

obtain as initial vector v 0 = 

2 

l 4 7 

+ i 

4 

116; 

= 

2 + 4/ 

4 + 1 6/^ 


#45: spirall ( [ 8 , -14, 7], t, [1 + i. 



- 0.9375 
- 1 . 875 


2 + 4 • i , 4 + 1 6 - i ] ) 

#46 : [2 A t, 4 A t] 

Notice, that the blue part of the parabola is the 
parametric curve corresponding to the Binet- 
representation. 

If we had chosen X 2 = 0.5, A = 0.25 , then the 
corresponding parametric curve with a suitable 
initial vector would be the arc of the parabola 
connecting z = 1 + / with z = 0 . 

To complete our tour through conic sections we 
will give now the construction of a 5/net- 
representation, which generates a parametric 
curve, which is part of a hyperbola. 


- 2.8125 


For a hyperbola we choose Ai l =l,A 2 >0, and A = — • Then p = 1 , q = -r , and 

A 

r = l + A 2 +A 3 = l + A 2 +— . The next figure is with A = f A = 2, and A = 0.5 , so 

A 
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p = 1, q = -3.5, r = 3, 5 . Since the parametric curve b(t) = 2‘ + i • 2 ' is part of the hyperbola 

as initial vector. 


given by v = — , we take v 0 = 

rn 


r i ^ 

2 

+ i • 

0.5 

X 

7J 


(0-25 y 


#47: spirall ( [ 1 , -3.5, 3.5], t, [1, 2, 4] + i - [ 1 , 0.5, 0.25]) 
#48: [2 A t, 0 . 5 A t] 


#49: pointseql ( [ 1 , -3.5, 3.5], [0, 6], [1, 2, 4] + i - [ 1 , 0.5, 0.25]) 
#50: We omit the list of points, which is the result of #49: 



So far, we have been able to construct special conics, and the question is: Is it possible to gen- 
erate any conic section as the interpolation of a /? //^/-representation of a suitable recurrence. 
The answer is affirmative, and the next two simple propositions are the key to the proof of the 
general result. 

Let / : C — > C be a function. We can apply / to the components of the initial vector v 0 to 


obtain the modified initial vector 
(0.14) 


f 



7(a)", 
/(a) • 
7(A)y 


If we use this modified initial vector in (0.3) we obtain a modified recurrence < z' n [ as re- 

( 1 «gN 

suit. In some particular cases the impact of / on the recurrence is the same as on the initial 
vector, i.e.: z' n = ffz n ) for all /ieN. 


Proposition 3.1 

Let | z' n > be a recurrence with initial vector (0.14), then 

l JnsN 

(1) f(z) = a- z,a eC, implies z' n = / (z n ) for all «eN 

(2) p + q + r = 1 and f(z) = a- z + b with a,b e C implies z' n = / ( z n ) for all «eN 

(3) p,q,r el and f(z) = a-z + b-z with aji e C implies z' n = f{z n ) for all «eN 

(4) p,q,r eM,, p + q + r = 1 and f(z) = a -z + b -z +c with a,b,c e C implies 

z' n = /(A) f° r al1 n e N 
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Proof: (1) is an evident consequence from the definition. To prove (2) observe 
(0.15) 



1 

0" 

' a- z n +b ^ 

f 

0 

0 

l 

a • z n+l + b 

= 

[p 

q 

q, 

K a-z n+1 +b) 

V 


a ' z «+i + b 
a • z n +2 + b 


p -(a ■ z n +b) + q - (a ■ z n+1 + b) + r fa ■ z n+2 +b)J 


a - z „+i + b 


f a ' Z n + l +b ^ 


(A z „ +1 f 

a-z n+2 +b 

= 

a ‘ Z n + 2 +b 

= 

f( Z n + 2 ) 

K afp-z n +q-z n+l +r-z n+1 ) + (p + q + r)-b)) 


a-(p- z n +q- z „ + 1 +r- z n+2 ) + b)j 


l/(U +3 )y 


The proof of (3) is by calculation 


P‘(a-z n + b-z n ) + qfa- z n+] +b- z n+l ) + rfa ■ z n+2 + b ■ z n+2 ) = 

(0-16) = afp-z„+q ■ z„ +1 + r ■ z n+2 ) + b ■ (p-z n +q- z B+1 + r ■ z n+2 ) = 

= a-z n+3 +b-z n+3 =/(z n+3 ) 

because conjugation of a complex number and multiplication with a real number commute. 

The proof of (4) is a combination of (0.15) and (0.16). This completes the proof. □ 

We changed only the initial vector, but the eigenvalues are kept, so the next proposition is a 
simple consequence of the last: 

Proposition 3.2 

Under the conditions of proposition 3.1 the Binet-representation of the recurrence 
with initial vector v 0 ' is b'(t ) = / (b(t)) . 

Proof: We have already proved, that b'(t) = f(b(t)) holds for all «eh. The coefficients in 

the /?me/-decomposition v () ' = a'e' + a 2 e 2 + a 2 e 2 \a,' e C 3 are unique, thus the proof is fin- 
ished. □ 

We summarize our results: 

• Given an arbitrary conic section T in C, we can construct a suitable recurrence and 
build a suitable initial vector v 0 so, that the resulting sequence {z„} )ieN is a discrete 
subset of T . 

• In the case of a circle or an ellipse we con construct recurrences with suitable initial 
vectors and {z„} neN c T , which are either periodic with a given period length in e M , 
or dense in T . 
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Proof: We gave explicit constructions for circles, a parabola, and a hyperbola. Since all cir- 
cles and all parabolas are similar, we can use proposition 3.1 (2) to map the special circle at 
any other circle or parabola by a similarity map z — > a ■ z + b with a,b eC. By proposition 
3.1 (4) a circle and a hyperbola can be mapped at any arbitrary ellipse or hyperbola by an af- 
fine map z — > a-z + b-z +c with a,b,c e C . This proves the first statement. 

To prove the second, we just need to start with a an eigenvalue of the fonn A = e' a n and 
a e Q , as in the example above, for a periodic orbit, and a e M \ Q for a dense orbit. □ 

This is the complete answer to Josefs question. 


(4) The Case of Multiple Eigenvalues 

Until now we have considered only examples with three distinct eigenvalues. It is even pos- 
sible to give a /fmet-representation in the case of multiple eigenvalues. In the following 
propositions M is always the matrix defined in (1.4). 

Proposition 4.1.1 

Let 2, be a simple and /t 2 be a double eigenvalue of the matrix M. 

[U fl] f O'. 

Set e 1 = Ay , e 2 = A 2 , and v = 1 . Then the following statements hold: 

j v^-2 y v^2 j 

(1) (e x ,e 2 ,vj is a basis of C 3 iff A x it A 2 . 

(2) ( M - A 2 ) v = e 2 , or equivalently: Mv = e 2 + A 2 v 

(3) M"v = n ■ Af 1 • e 2 + Af ■ v for n > 0 


Proof: (1) An easy calculation gives det(ej,e,,v) = (A 2 -A l ) 2 *0oA 2 *A l . 

(2) We expand (z - A l )(z - Af) 2 = z 3 -(A, +2Af)-z 2 +(2A l A 1 +A 2 2 )-z-A l Af and see, that the 
parameters in the matrix M are p = A l Af , q = -(2 A l A 2 + Af ),r = A l + 2 A 2 . We proceed: 


(~*i 

(M - Af)v - 0 

V P 


1 

-T 2 

q 



~A 2 + 2 A 2 


— (2/Ij A 2 + A 2 ) + (Tj + Af) • 2 A 2 



(3) As a consequence of (2) we have: M 2 v = M(e 2 + A 2 v) = A 2 e 2 + A 2 Mv = 2 A 2 e 2 + Afv , and by 
induction M"v = n ■ Af • e 2 + Af ■ v . □ 


We notice, that this proposition is true, even if one eigenvalue is zero. The next one settles the 
case of a triple eigenvalue. 

Proposition 4.1.2 

m r o a roA 

Let A it 0 be a triple eigenvalue of M and e= A , v, = 1 , and v 2 = 0 . Then: 

, A 2 , 2 A , T, 
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(1) (e,v x ,vf) is a basis of C 3 . 

(2) ( M - A)v x = e or equivalently: Mv l = A ■ v x + e 

(3) ( M - a)v 2 = Vj or equivalently : Mv 2 = A ■ v 2 + v l 

Proof: (1) We calculate det(e,v 1? v 2 ) = 1^0. 

(2) Since A, is a triple eigenvalue, we find for the parameters in M: p = A 3 ,q = -3 A 2 , r = 3 A . 

( -A 1 0 V 0 W 1 W 1 ^ 

So ( M-A)v= 0 - A 1 1= -A + 2A = A =e follows. 

y A 3 -3 A 2 2A j[2A j {-3A 2 +4A 2 J [a 2 , 

(3) A similar calculation proves (3). □ 

Now we have all pieces together to state 

Proposition 4.1.3 

Let the notations be as above and let w := a x ■ e x + a 2 ■ e 2 + a 3 ■ v, a k eC for k = 1, 2, 3 . 
Then in the following formulae are valid: 

In the case of a double eigenvalue A 2 : 


(1) M n w = a x Af ■ e x + (a 2 • Af + a 3 ■ n ■ Af 1 ) • e 2 + a 3 ■ Af ■ v for n e N c 


(2) If A 2 ^ 0 , then M"w = a x Af ■ e x + Af |^a 2 + a 3 ■ — J -e 2 +a 3 - v for «eH 0 

a x ■ e x + a 2 ■ e 2 + a 3 ■ v for n = 0 

(3) If A 2 =0 , then M n w-< a 1 A l -e 1 +a 3 -v for n = 1 

a x Af ■ e x for n> 2 

In the case of a triple eigenvalue A^O : 


(4) M"v j = A" ■ Vj + n -A" 1 • e for n > 0 

(5) M"v 2 = A" ^ ” < 2 -1) • e + n ■ v x + v 2 j for n > 0 

(6) M n w = A n [{na x + ”f 1] a 2 + a 3 ) ■ e + (a x + naf)- v x +a 2 -v 2 ^j for n> 0 


Proof: (1) is a consequence of proposition 4.1.1 part (3). (2) is equivalent to (1) if A 2 ^ 0 . 
(3) If A 2 = 0 , then 


o l oYi 


0 1 ovo 


Me 2 =00 1 0 = 0 , Mv = 0 0 1 0 = 1 , and M 2 v = A x =e x . 


0 0 A Jlo 


0 0 A1 


Calculating Mw and M 2 w proves (3). 

(4) From 4.1.2 part (2) we have M 2 Vj = A 2 +A(A-e + A-v x ) = 2 A 2 ■ e + A 2 ■ v x and (4) follows 


by induction. (5) follows from proposition 4.1.2 part (3) and the previous statement. 
(6) is the result of the linearity of the powers of M. □ 
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We choose the initial vector w = z x and pass in the previous proposition to the components 


of the resulting vector M"w . This leads us to 


Proposition 4.2 (Um^f-representation for multiple eigenvalues) 

Let A 2 be a double eigenvalue and let a l ,a 2 ,a : , be the coefficients of the Binet- 
decomposition of the in itial vector . 


(1) If A 2 0 , then z n = af" + a 2 Xf 1 + — J for n > 0 

(2) If A 2 = 0 , then z n = aff for n > 2 
Let A be a triple eigenvalue, A ^ 0 . Then 


(3) z n = a, ■ A" +n-a 2 ■ A "- 1 + ’AlLJl . ^ . r ~ 2 for n> 0 

Remarks: (1) We note, that only (1) and (3) admit a continuous interpolation in the sense, 
that we get z 0 ,z l ,z 2 for t = 0,t = l,t = 2 if we replace the discrete variable n by the continu- 
ous variable t e [0,oo) . 

(2) If we had chosen A = 0 as a triple eigenvalue, then we had z n = 0 for all n > 3 , whatever 
the initial values had been. 

(3) In the first case, convergence or divergence of the sequence can depend on the 

eigenvalues and on the initial vector. In the other cases (2) and (3) the convergence depends 
only on the eigenvalues f resp. A . 

(4) It is no problem, to write a DERLVE- function, which calculates the appropriate Binet- 
representations for all different cases. We present the DERLVE- code and two examples: 


auxl (p, t, v) ■■= 

Prog 

e >= EIGENVALUES ([ [0, 1, 0], [0, 0, 1], p] ) 
w ~ VECTOR (SUBST ( [1; r; r A 2], r, e;k), k, 1, 3) 

(SOLUTIONS ( [a, (5, y] -w = [vjl; vj2; vj3], [a, P, yl ) -VECTOR ( IF ( k # 0, k A t, 0, 0), k, e))jl 


aux2a(p, n, v) ■— 

Prog 

w := [1, 1, 0; pjl. Pi 2, 1; pil A 2, pi2 A 2, 2-pi2] 

(SOLUTIONS (w- [a, b, c] = v, [a, b, c] ) ) i 1 • [IF (p i 1 = 0, 0, pil A n), pi2 A n, n-pi2 A (n - 1)] 

aux2b(p, n, v) 

Prog 

w ■= [1, 1, 0; p;l, 0, 1; p i 1 A 2 , 0, 0] 

(SOLUTIONS (w- [a, b, c] = v, [a, b, c] ) ) i 1 • [p; l A n, 0, 0] 

aux2 (p, n, v) : = 

If p i 2 = 0 

aux2b(p, n, v) 
aux2a(p, n, v) 

aux3 (p, n, v) >= 

Prog 

w - [1, 0, 0; p;l, 1, 0; p;l A 2, 2-pU, 1] 

(SOLUTIONS (w- [a, b, c] = v, [a, b, c] ) ) i 1 • [p; l A n, n-pil A (n - 1), n-(n - l)/2-pil A (n - 2)] 
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binettp, n, v) == 

Prog 

s ■= SOLUTIONS (x A 3 = p-[l, x, x A 2], x) 
w ■= I F ( s | 1 + 2 • s i 2 = p I 3 , S, REVERSE (s)) 
1 = DIMENSION (s) 

If 1 = 3 

auxl (p, n, v) 

If 1 = 2 

aux2 (w, n, v) 
aux3 (s, n, v) 


comptoreal ( z ) ■= [RE(z), IM(z)] 

pointseq(p, n, m, v) •= VECTOR (comptoreal (binet (p, k, v) ) , k, n, m) 


First example: 

We take the single eigenvalue \ = 0.6 + 0.8/ and the double eigenvalue \ = 0.6 - 0.8/ . This 
gives p = 0.6 - 0.8/, q = -\.12 + 0.96/, r = 1 .8 - 0.8/ . We omit the output of the resulting list of 
points but print only the command and the graphical representation of the result: 


pointseq ( [ 0 . 6 - 0.8-i, -1.72 + 0.96-i, 1.8 - 0.8-i], 0, 10, [1, i, 1 + i]) 



We notice, that the sequence of points is divergent. 
The reason is, that the factor 


f \ 

1 + A 


V 


JL 


in 


2 y 


Z ;i — 

tends to infinity, though |/tj 


V ^2 J 


is unbounded, if n 
"| = U 2 "| = 1 for all n . 


Second example: 

This example works with \ = 0.5 + 0.8/ and 
\ = 0.5 - 0.8/ , so the absolute values are just be- 
low 1: 


pointseq ( [0 . 445 - 0.712-i, -1.39 + 0.8-i, 1.5 - 



-A 


0.8-i] , 0, 100, [1, i, 1 + i] ) 


We can see, that the values are at first 
increasing, but since 

r \ 


lim a x / lj" + a 2 A 2 


1 + - 


A 


= 0 , the se- 


'2 J 


quence of points is convergent to the 
origin. 


Finally we look at limits of quotient 
sequences l > , because the Bi- 

1 n J hgN 

«e/-representations allow us to give 
exact conditions on the convergence. 
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As before, we must distinguish three cases: 


(0 


(ii) 


(iii) 


The characteristic equation of the recurrence has three distinct roots. Then, 
roughly speaking, the eigenvalue with largest absolute value, which is present in 
the Binet-representation, wins: Let \f\ > |/l 2 | > \l 3 \ , then by proposition 2.2 




n + 1 


+ cc 2 A 2 


n + 1 


+ ( 2 3 A 3 


n + 1 


a x Af + af if + a 3 Af 


Since 


A 

< 1 and 

£i 





< 1 , the limit is 


lim z,,„ ^ 

n ^°° A l n (a l + a 2 (A 2 " /A l ) + a 3 (A 3 /A 1 )" s J ' 


The characteristic equation has two distinct roots. Then again the eigenvalue with 
largest absolute value, which is present in the Binet-representation, wins. If the ei- 
genvalues are different in absolute value, then let A u be the eigenvalue with lar- 
ger absolute value. We find by proposition 4.2 and similar considerations as 


<afj Aj" +I + a 2 A 2 


n + 1 


above, that lim 


1 + 


n + 1 


— = lim - 

7 n — >oo 


A. 


2 ) _ 


= A„ 


a x Af + a 2 Af 


1 + 


V 


A. 


2 J 


In the case of a triple eigenvalue we get by proposition 4.2, that 

1 n +l i in, n{n+\) ^ i n - 1 


, . z„ , , , . a. - A 

11111 -^ = Inn 1 


+ n-a 2 - A" + : 2 


■a 3 -A n 


»->« a x ■ A" + n ■ cc 2 ■ A n ~ l + ^-a 3 - A"- 2 


„ , . cn, ■ A + n • cx n • A + 

= A ■ lim 


n(n+ 1) 
2 


a , 


n(n- 1) 


n ^ co a l - A +n-a 2 -A + 2 


a , 


Comments: In the case of two or three distinct eigenvalues a limit must not exist, while in 
the case of a triple eigenvalue a limit always exists, if A ^ 0 . If 2 = 0, then z n = 0 for n > 3 

and the quotients are not defined. Note, that the presence of a particular eigenvalue in the Bi- 
net-representation depends on the initial values. 

We finish with some number theoretic considerations: The ffinet-decomposition of the initial 
vector and the ZimeZ-rcpresentation together show, that the following proposition holds: 

Proposition 4.3 

Let F(p,q,r,z 0 ,z l ,z 2 ) be an intermediary number field Q c: F(p,q,r,z 0 ,z l ,z 2 ) a C, 
which contains the coefficients p, q, r of the recurrence and the initial values z 0 ,z,,z 2 . 
Then the following statements hold: 

(V { Z n} neN ^ F (P’<l’ r ,Zo’ Z l’ Z 2) 

£ 

(2) If one of the limits lim z n or lim— ^ exists, then this limit belongs to the unique 

X— >00 X — >00 7 

field extension F of F(p,q,r,z 0 ,z l ,z 2 ) of at most degree three, which is generated 
by the roots of the characteristic equation of the recurrence. 

Proof: (1) The definition of the recurrence uses only the four rational operations, so all the 
z n stay inside F{p,q,r,z 0 ,z l ,z 2 ) . (2) The eigenvalues are roots of the characteristic polyno- 
mial, so a standard argument from Galois- theory says, that they belong to an extension field 
of at most degree three. The coefficients in the B in et-decomp o sitio n are solutions of a system 
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of linear equations in this extension field, so they must belong to this field too. The limits, if 
they exist, depend only on the coefficients of the 5me/-decomposition and on the roots of the 
characteristic equation. So they belong to the extension field. □ 

Let us consider Carl’s and Josefs example once again: The initial values are pairs of rational 
numbers, or in complex notation, numbers with rational real and imaginary parts. The same is 
true for the eigenvalues (2.7). As a consequence, the limit z K = 0.4 + 0.4/ has rational real and 
imaginary part. 

(5) A Theorem on Coefficients 


We only mention, that everything presented here can easily and in a straightforward manner 
be generalized to higher dimensional recurrences 

m 

Z „ +m = Yj Q k Z n + k With illitial ValueS Z 0> Z l>-> Z m - 1 G C ' 


Thus we state the next propositions in full generality. The convergence of the recurrences 
depends generally on the eigenvalues, they depend on the coefficients a v a 2 ,...,a m of pa- 
rameters and on the initial vector v 0 = (z 0 ,z 1 ,...,z m ) . In a very particular case the convergence 
is guaranteed by a condition on the coefficients of the recursion. To state the proposition pre- 

f k if a. + 0 

• „ _ i _ . • x a _ _r„ii : cv \ . J * 


cisely, we introduce the following notation: 8{a k ) := 


Proposition 5.1 


0 if a k = 0 


n— i 

Let p(z ) :=^a k z k be a polynomial with nonnegative coefficients a k > 0 for 

k = 0 

n — 1 

k = 1 -l, and ^a k =:s . Then the following statements are true: 

k = o 

(1) If s = 1 and z is a solution of the equation z" = p(z) , then z| < 1 . In particular 
Zj = 1 is a solution. 

(2) There exists a positive solution q of the equation z” = p(z) , and for all other so- 
lutions z we have |z| < . 

(3) If s < 1, then there exists a positive solution c < 1 of the equation z" = p(z) , and 
for all other solutions z we have Izl < £ < 1 . 


(4) If all coefficients are nonnegative, 5 = 1, and m := gcd(£(a 0 ), S(a 1 ),..., S(a n _f)) , 
then the equation z" = p{z) has exactly m solutions z.,j = l,...,m,with | z ; . | = 1 . 
All other solutions satisfy Izl < 1 . 


Proof: (1) We assume, that there is a solution z with z > 1 . We obtain the following chain 

n — 1 n - 1 

of equalities/inequalities: \z n | = ^ a k z k < ^ \a k | • |z A ' | < 5 • |z H_1 1 = |z" _1 1 <=> |z| < 1 . This is a 

k = 0 k = 0 

contradiction to |z| > 1 and the first part is proved. Clearly z = 1 is a solution because we have 

5 = 1. 
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(2) We first show, that exactly one positive solution exists. Consider g(x) := — . After 


possibly cancelling equal powers of x we have: g(0) = 0 . On the other hand limg(x) = oo , 

x — >oo 

so, by continuity, there must be a real number £ > 0 with g(£) = 0 . Now let z be another 

( z V n- 1 a f z \ k z 

solution, then division by <£" gives — = V — k -r — . We substitute 6 := — and obtain: 

J 77 ° p pn-k p p 

J k=o b v S y S 

6*" = ■ <9 = 1 is a solution of this last equation, so by proposition 1.1: ^ = 1 ■ 

k=0 4 k = 0 £ 

2 I i 

Since d, > 0 , we find by ( 1), that — < 1 or equivalently: z < £ . This also proves the unique- 
ness of £ . 

(3) We only need to prove, that £ < 1 if 5 < 1 . Assume, that £ > 1 . Then \ , and as a 

consequence 1 = V — ^ <V — <— and equivalently £ < s < 1 . This ia a contradiction, so 
*=o <T“ *=o I £ 

we have proved < 1 . 

(4) If the gcd equals 1, then the statement is a theorem of Ostrovsky [2, p. 3], which can be 
proved along the lines of (1) and (2). The condition on the the gcd has the following interpre- 
tation: Let m := gcd(S(a 0 ),S(a l ),..., 8{a n _ l )) , then z” - p(z) = q(z m ) is a polynomial of z m 

L>-1 L>~ 1 

and q(z) := z v - ^ b k ■ z k . Obviously the coefficients satisfy b k = a m . k and ^ h k = 1 . Now we 

k = 0 k = 0 

have gcd(^(/i 0 ),..., 8{a u _ x ),v) = 1 , so by Ostrovsky’s theorem z = 1 is the unique solution of 
q(z) = 0 with |z| = 1 . Thus the numbers z . := e 2m ]/m ,j e {0,1, ...,/« -1} are the only solutions 
of z” = p(z) with |z| = land the proposition is completely proved.n 

The following corollaries are immediate consequences of proposition 5.1. 

Corollary 5.2 

m 

Let z n+m = Z«A- be a recurrence with initial values z 0 ,z x ,...,z m _ x eC, nonnegative 

k= 1 

n — 1 

coefficients a k , ’ s ^_ i a k =1, arid gcd(8{a () ),8{a l ),...,d(a n _ x )) = 1 . Then the sequence 


iA L„ converges. 


Proof: The conditions stated above guarantee, that there is exactly one simple eigenvalue 
\ = 1 and all other eigenvalues satisfy \\ | < 1, k e {2,3,...,n} . So even in the case of multiple 

eigenvalues the corresponding fh'net-representation is convergent, and the convergence is 
independent of the initial values, because the powers of the eigenvalues tend to zero as the 
exponents tend to infinity, for k e {2,3,...,«} , and to I for \ = 1 .□ 
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Corollary 5.3 


Let z n + m be a recurrence with initial values z 0 ,z l ,...,z m _ l eC, nonnegative 

k = 1 

n — 1 

coefficients a k , and 'f^a k = s < 1 . Then the sequence {z n \ nrU converges to zero. 

k = 0 

Proof: By proposition 5.1. (3) all eigenvalues satisfy \A k \ <1, and consequently all powers 
in the corresponding fime/-reprcscntation converge to zero whatever the initial values are. □ 

Remarks: (1) This corollary was proved for n = 3 with different means as proposition 1.2 . 
Josefs example from [1, p. 5] with p = a 0 = q = a x = 0.5 and r = a 2 = 0 satisfies the condi- 
tions of corollary 5.2: p + q + r = 1 and gcd(0, 1, 3) = 1 , so the sequence is convergent with 
every initial triangle. 

(3) The condition on the ged in corollary 5.2 cannot be omitted, because this condition on the 
ged is necessary to guarantee convergence, as the following example will show. Let n = 3 , 
a l = 1 , and a 2 =a 3 = 0 . This is the case of the triple eigenvalue A = 1 . We have the Binet- 
representation: 


i« . ih- 1 i 1) irt-2 

z n = a l ■ 1 + n ■ a 2 ■ 1 -i a 3 ■ 1 


, n —2 n(n- 1) 

a 3 ■ 1 = a l +n-a 2 - 1 a 3 


So lim z n = oo if a 2 and a 3 do not vanish both. It is possible to construct examples for n = 4 

n—> oo 

which are bounded or even periodic but not convergent. 

(6) Conclusions 


A small shift from points in the plane to complex numbers allows us to use the machinery of 
linear algebra. Especially the use of eigenvalues of the matrix, which describes the recurrence, 
allows us to understand some geometrically surprising phenomena occurring with computer 
experiments, such as points on a circle or ellipse. The corresponding Binet - representation al- 
lows to build something like a construction kit for recurrences with a prescribed behaviour. 
The /?me/-represcntation gives us complete information about convergence or divergence of 
the sequence of points defined by recurrences. The convergence/divergence is governed by 
the eigenvalues, resp. the roots of the characteristic equation and all we have to look for, is 
lim R(n) ■ A " , where R(n) is a rational function of n . Finally a CAS, like DERIVE, allows us 

n—> co 

to visualize the theory of recurrences. 
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Our long time DUG Member Robert Setifsent a mail including a picture of the famous Man- 
delbrot Set - the Apple Man: 

Bonjour. 

Can Derive program a function that creates an image such as the attached file? 

Thank you and best regards 



Robert 

It is mere chance that I had also an exchange of mails with Rob Gough concerning the 
famous Mandelbrot set. He sent a DERIVE file and complained about boring calculation 
times. 


Dear Josef, 

However, of greater importance to me is my ability to program in Derive (for my own pleas- 
ure). With a lot of difficulty (because the Help functions don’t work properly) I managed to 
create a program to produce the Henon Attractor (Chaos Theory). So my next task was to 
produce the Mandelbrot set. This followed on relatively easily from the previous work, but I 
could not produce a graphical presentation of this. I know how to produce the data but can- 
not display it! I produced a 2D iteration based on x:=x A 2 -y A 2+a, y:=2*x*y+b that produced a 
3D vector [a,b,c] where c is a measure of the stability of the iteration. My intention was to use 
the third axis to represent the colour of the 2D [a,b] plot. Can I do this? If it would help I could 
enclose a pdf file of the work I have done. 

Best wishes 


Rob 
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The Mandelbrot Set 

Rob Gough 


The Mandelbrot Set defines a region where a certain recursive function remains stable and 
does not produce values heading for infinity. An outline is shown in 'Chaos Under Control' by 
D Peak and M Frame, 1994. The shape is mainly confined to a region with -2<a<+0.5 on the 
horizontal axis (a-axis)and -1<b<+1 on the vertical axis (b-axis). The aim is to concentrate on 
any part of this 2D plot within the limits [A1 , A2, B1 , B2] where A1 <a<A2 and B1 >b>B2 with 
any desired pixel detail ( ). To create a linear sequence of coordinates suitable for the De- 
rive VECTOR function it is necessary to scan the desired region in the fashion of a CRT dis- 
play from top left [A1 ,B1] to bottom right [A2,B2], 


A2 - A1 B1 - B2 

K := 1 + , L := 1 + , 3 := K-L 

8 5 

[k :e Integer, 1 :e Integer, j :e Integer] 

[0 < k < K, 0 < 1 < L, 0 < j < 3 - 1] 


1 := FLOOR 


k := j — 1 ■ K 


a::Al+ S- j - FLOOR 


B1 - 8. FLOOR 


#6: [a := A1 + S-k, b := B1 - 5-1] 

The 2-dimensional [a,b] plot can be coloured with reference to the stability criterion of the 
recursion. This is based on five numbers [Cl, C2, C3, C4, C5] where C1<C2<C3<C4<C5 
with Cl the least stable and C5 the most stable. For example, if the recursion is still within 
certain limits after C5 iterations then the parameters [a,b] is considered stable and the third 
dimension records '5' and the coordinates [a,b,5] are returned. 

The recursion is based on the dummy variables [x,y]: 

x:= x A 2 -y A 2 +a 

y:= 2*x*y + b 

and the stability criterion is based on x A 2 + y A 2 <4. 

#7: [x :=, y :=] 

XY_AUX(A1, A2 , Bl, B2, 8, n, x, y) := 

If n = 0 
#8: [x , y] 

XY_AUX(A1 , A2 , Bl, B2, 8, n - 1, x A 2 - y A 2 + a, 2-x-y + b) 

#9: XY(A1 , A2 , Bl, B2, 8, n) := XY_AUX(AL, A2, Bl, B2, 8, n, 0, 0) 


#10: DMAX(A1 , A2, Bl, B2, 8, n) := XY(A1, A2, Bl, B2, 8, n) 
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#11: levsQ := [2, 5, 10, 20, BO] 

M(M, A2 , Bl, B2 , 5, Cl, C2, C3, C4, CB , a, b) ;= 

Prog 
n 0 
Loop 

If 0 < n < Cl a DMAXCAl , A2, Bl, B2, S, Cl, C2, C3, C4, C5 , a, b, n) > 4 
RETURN [a, b, 0] 

If Cl < n < C2 a CMAXCA1, A2, Bl , B2, 6, Cl, C2, C3, C4, CS , a, b, n) > 4 
RETURN [a, b, 1] 

#12: If C2 < n < C3 a DMAX£A1, A2, Bl, B2, 5, Cl, C2, C3, C4, C5 , a, b, n) > 4 

RETURN [a, b, 2] 

If C3 < n < C4 a DHAX(A1, A2, Bl, B2, S, Cl, C2, C3, C4, CS , a, b, n) > 4 
RETURN [a, b, 3] 

If C4 < n < CS a DMAXfAl, A2, Bl, B2, 5, Cl, C2, C3, C4, CB , a, b, n) > 4 
RETURN [a, b, 4] 

If n = CB 

RETURN [a, b, B] 

n := n + 1 


#13: MB(A1 , A2, Bl, B2, 8, Cl, C2, CB , C4, C5 , j) :=M(A1, A2, Bl , B2, 8, Cl, C2 , C3, C4, C5 , a, b) 

#14: MBS(A1, A2, Bl, B2, 8, Cl, C2, C3 , C4 , C5) := VECTOR(MB(Al , A2, Bl, B2, 8, Cl, C2, CB , C4, CS , j) , j, 0, J - 1) 

Lastly, these 3-dimensional coordinates [a,b,c] are returned as a vector sequence of length 
J. 

Here is a more detailed example: 


#25: 

[A1 := -1.3, 

Bl := 

0.4, 

A2 := 

-0 

■ 7, 

B2 := 

-0.4, 

5 

:= 0.02] 


#26: 

[Cl := 5, C2 : 

1 1 

K 1 

o 

C3 

:= 20, 

C4 

:= 

LTi 

U 

O 

LTl 

:= 100, 

i 1 

O 

LTl 
J [ 

c 


#27: 

MBS (-1.3, -0 

.7, 0 

■ 4, 

-0.4, 

O.i 

02, 

5, 10 

, 20, 

50 

, 100) 


#28: 

i 1 

i 1 

1 

K 1 

UJ 

o 

1], 

[-1. 

28, 0 

■ 4, 

1] 

, [-1. 

26, 0 

■ 4, 

2], [-1.24, 0.4, 

2], [-1.2 


K 1 

i i 

i — i 
1 

K 1 

K 1 

r-j 

, 0.4 

, 1] 

, [-1. 

.1, 

0. 

1 1 

1 — 1 

’tf- 

[-1. 

08, 

0 
1 

1 

1 1 

1 — 1 
1 — 1 

o 

i 1 

1 — 1 

T T 

O 


#28 consists of 1271 triples. 

My first answer was as follows: 

I introduce variable setl for the list of points (which is not printed here). 

#26: setl := MBS (-1.3, -0.7, 0.4, -0.4, 0.02, 5, 10, 20, 50, 100) 

#27: setl := [[-1.3, 0.4, 1], [-1.2S, 0.4, 1], [-1.26, 0.4, 2], [-1.24, 0.' 

[-1.14, 0.4, 1], [-1.12, 0,4, 1], [-1,1, 0.4, 1], [-1.0S, 0,4, 1], 

Calculation of #26 needs a very long time (7000 seconds! !). I believe that the reason is the recursive 
nature of # 8. 1 am quite sure that this should (and could) be improved. 

This is the trick: The SELECT-command separates the points according to the depth of iteration and 
extracts the first two columns (x,y) making the points ready for plot. The resulting matrices are plotted 
in different colours. Make sure that the points are plotted discrete and you may change the point size. 

setll := VECTOR((SELECT(v = k_, v, setl))u[l, 2], k_, 0, 5) 

3 


#2S: 
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Expression #28 can be plotted in the 2D Plot Window. It is not necessary approximating #17 but take 
care that in the Plot Window Options > Approximating Before Plotting is activated. Each Plot results 
in another combination of colours. 



Now, the 3D-presentation: I apply SELECT again, but the result does not lead to a satisfying plot 
(#29), because all points are connected. 

setlll := VECTOR (SELECT (v = k_, v, setl) , k_, 0, 5) 

#29: 3 

setlSd := VECTOR (VECT0R([v], v, setlll), k, DIM (setlll)) 

#30: k 

The command in #30 produces single points collected according to their levels (0 - 5) in matrices 
(lists) which are ready for plotting in 3D. 



Point Size Large (Left) and the Top View with Point Size Medium (right) 


As the readers might assume I was not really satisfied with my first answers. Not because of the result- 
ing plots but because of the enormous calculation times. I tried to do better which was not really nasty 
for me because I have been interested and fascinated by all kinds of fractals (see earlier Newsletters). 

I was sure that I had to overcome the recursive procedure by writing a program including all necessary 
parameters in the parameter list of the function (program) to be developed. 
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I define some “levels of escape”. pos() relates the number of the “escaping iteration” to the position 
of the escape level in the list. Given 1 evsl and iteration #45 escapes, then pos(l evsl , 45) = 3 and 
pos(levs5 ,45) =4 but pos(l evs2 , 45) = 1, etc. 


#1: 

levs : 

= [2, 5 

, 10, 20, 

50] 


#2: 

levsl 

:= [5, 

10, 20, 50 

100] 

#3: 

levs 2 

0 
l n 

1 i 

1 1 

100, 130, 

170, 

200] 

#4: 

levs 3 

:= [100 

, 130, 170 

200 

300] 

#5: 

levs4 

:= [5, 

10, 20, 50 

75, 

I 1 

□ 
LTi 
1 — 1 

o 
o 
1 — 1 

#6: 

levs 5 

:= [5, 

10, 20, 40 

50, 

60, 70, 80, 90, 100] 

#7: 

pos(l. 

n) := 

POSITION (true, 

VECTOR (m - n > 0, m, 


cri t() is the core part of my procedure. It returns the initial point [a,b] together with its “escaping 
number” - depending on the level list 1_ - which is later associated with a plot colour. 

crit(a, b, 1_, zO, zl, n, n_) := 

Prog 

zO a + l ■ b 
n := MAX (levs) 
n_ := 1 
Loop 

#7 : If n_ > n 

RETURN [a, b, DIM(1_)] 
zl := zO a 2 + a + b-i. 

If ABS(zl) > 2 

RETURN [a, b, pos(l_, n_) ] 
zO := zl 
n_ : + 1 

Expression #8 scans the region for al < x < a2 and bl <y < b2 in S-steps. 

#8: MS(al, bl, a2, b2, 8, 1_, zO, zl, n, n_, x_, y_) := 

APPEND(VECTOR (VECTOR (crit (a, b, 1_, zO, zl, n, n_) , a, al, a2, 8), b, b2, bl, 8)) 
Let’s see if it works and if there is really a success in speeding up the calculation: 

#11: set := MS(-1.3, 0.4, -0.7, -0.4, 0.02, [5, 10, 20, 50, 100]) 

needs only 43 seconds!! 


I don’t print the set set . In fact I removed it for saving memory! 

VECTOR ( (SELECT (v = k_, v, set))u[l, 2], k_, 5) 
#12: 3 

The nice thing is zooming in. So I zoom into the black rectangle top 
right. 

To make it easier I invented a function for plotting. 
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plot(s_, 1_) := VECTOR((SELECT(v = k_, v, s_))u[l, 2], k_, DIM(1_) + 1) 
#13: 3 


#14: SEHHHBEES1 

This gives the same picture as above. Then I deleted the point list set in order to save memory. 


#15: set := 

What is following was all done using this procedure. But at a certain moment I considered that 
formatting the out put and performing the output needs also calculation time and memory. In fact I do 
not need the list of the points, I only want them to be plotted. So I wrote a very short program 
combining calculation and plotting as well: 

MprogCal, bl, a2, b2 , 5, 1_, zO, zl, n_, n, x_, y_, set) := 

Prog 

#16: set := AF'F J END(VECTOR (VECTOR (cnt ('a, b, 1_, zO, zl, n_, n) , a, al, a2, 8), b, b2, bl, 8)) 

VECTOR C (SELECT CVi 8 = k_, v, set))u[l, 2], k_, + 1) 


Using Mprog I started the journey into the interior of the boundary region: 


#17: Mprog (-1 . 2 , 0.35, -1, 0.2, 0.005, levsS) 

needs 44 sec fior calculating and plotting 



#1S: 


Mprog (-1 . 13 , 0.3, -1.1, 0.2S, 0.0005, levs5) 



#19: Mprog (-1.127, 0.287, -1.121, 0.28, 0.00005, levsS) 



Same detail plotted using different colour palettes. 
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Performing the last zoom 

#20: MprogC’-l . 123 , 0.2S2, -1.122, 0.281, 0.00001, levsS) 



And tried this function for plotting the whole “Appel Man”. . . 

. . . followed by a short visit of the famous “Sea Horse Valley” 
#21: MprogC-2, 1.2, 1, -1.2, 0.02, levs5) 



#22: MprogC-0 . S57S, 0.1692, -0.6623, 0.0225, 0.005, levs5) 
#23: MprogC-0. S57S, 0.1692, -0.6623, 0.0225, 0.0025, levsS) 








p 


30 
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#24: MprogC-0. 75075, 0.10816, -0.742, 0.09898, 0.0001, [5, 10, 20, 30, 40, 50, 60, 70, 80, 90, 100, 150]) 



#24: Mprog(-0. 747, 0.1057, -0.745, 0.1046, 0.00001, levs5) 



Did you ever try changing the iteration zl:=zO A 2 + a + b-i to zl:=zO A 4 + a + b-i or zl:=zO A 6 + a + b-i 
or zl :=zO A 4 + zO A 2 + a + b-i or . . . 



See more spectacular Mandelbrot sets on page 50. 
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Comments on Klaus Korner’s Contribution in DNL#83 

David Halprin, Australia 

Josef 

Hello from downunder. 

Even though I haven't written to you for a long time, my interest in DUG had never 
waned. 

Your last DNL #83 grabbed my extra-special attention with the correspondence to and 
fro with Klaus Korner. My interest in the aerofoil profile started about 35 years ago when I 
came across a book by Richard von Mises in the university library. At the time there were no 
personal computers and I wondered how to plot such a curve. With the advent of the IBM PC 
locally in 1987 I still had trouble working out the right approach. 

However, 'the penny dropped' when I was working with curves, that were de- 
rived/transformed from other curves, such as:- Bertrand, Caustics (Dia-, Cata-), Caylean, 
Centroid (Centre of Gravity) Lines, Cissoids, Conchoids, Cycloids (Epi-, Hypo-, Prolate-, Cur- 
tate-), Evolute & Involute, Glissettes, Reflective Inverse, Isoptics, Orthoptics, Mannheim, Or- 
thotomic, Orthogonal Trajectories, Pedal, Contra-Pedal, Paragraph (Parallel, Telegraph etc.), 
Radial, Pursuit, Negative Pursuit, Pseudo-Pursuit, and many other pseudo curves u.s.w.. 

Here was my initial interest in going backwards and forwards from Involute to Evolute 
and all the other such pairs of curves with their inverses, i.e. Base Curve to Derived Curve, 
then Derived Curve back to Base Curve. 

Hence, the "AHA Moment" had arrived for the Joukowski Aerofoil profile, it being yet 
another derived/transformed curve, and in that particular case, the base curve is the nomi- 
nated circle. 


_ , . 1 x - ly 

1 J w 2 2 

E t * * y 


Lat Z 2 - u + xv 



JOUKOWSKI PROFILE DEFINITION 
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Here are some graphs of the airofoils produced with JOUK.MTH and JOUKOW-1 .MTH. 


It was very interesting for me (Josef) that plotting with DERIVE XM was much easier than 
plotting with DERIVE 6.1 . See the gap in the 6.1 . aerofoils (first graph) and compare with the 
respective DERIVE XM-graph. Both were plotted with the same interval for the parameter in 
JOUK.MTH. 




The next plots are results from JOUKOW-1 .MTH. 


f 

( L,t 
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+1) + h ) + 1 ^ 

j 
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L ■ t 2 2 

h-(l + i) + e 'i/CCh + 1) + h ) - 1 


1 - 


i * t 2 2 

v h-(l + i) + e --/((h + 1) + h ) + 1 j 


IM 


i-t 2 2 

h-(l + i) + e 'i/CCh + 1) + h ) - 1 


1 - 


i-t 2 2 

v h- (1 + l) + e 'i/CCh + 1) + h ) + 1 j 


h = 1/8, 1/4, 1/2, 1, 3/2 and 2. Plotting with DERIVE 6.1 needs a lot of calculation time (first 
figure) 
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Plotting with DERIVE XM is a work of an instant! 
I don’t know what makes the difference? 



P.S. 28 October 2011 

I wrote this letter/paper yesterday afternoon, and before I went to bed I discovered an online 
eBook website that allowed me to purchase the 1959 Reprint of 'The Theory of Flight', so I 
shall be able to relive the experience after 35 years and see what von Mises' mathematics 
does. 

I shall send you a copy later. 

This is the link where David bought the copy of Mises’ Theory of flight: 

http://isbnbookfinder.eom/Theorv-of-Flight/p 179525/ 

Regards 
David Halprin 


These are some links which I found wrt David’s paper: 

http://math.fullerton.edu/mathews/c2003/JoukowskiTransBib/Links/JoukowskiTransBib Ink 3.html 
http://de.wikipedia.org/wiki/Konforme Abbildung 
http://www.vaxman.de/analog computing/joukowski/ioukowski.html 
http://www.luftpiraten.de/glos jOO.html 


David sent another paper which I will present in the next DNL. 

Many thanks for your persisting interest in cooperating with the DUG, Josef 
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Some requests concerning the TIs: 

From Hannes Scherzer, Austria 

I am owner of Tl-Nspire™ CAS Version 3.0.1 .1753. My question is: 

How to shade an area between two curves. The TI-83 PLUS offers the function shade(. . . .). 
In DERIVE I can enter: 0 < y < f(x) AND (x > 0 AND x < 4). 

Is there a possibility to do this with Nsipre, too? 

DNL: 

I don’t know how to do this. Is there any Nspire user who can help? 


From Nils Hahnfeld, Virgin Islands 

(1) How can I convince the TI89 that (-1) A n * (-1) A n = 1? 
Why does it refuse simplifying this? 

(2) Do you know a trick how to simplify (2n)!/n! or (n+3)!/n! ? 


DNL: 


DERIVE does not have any problems with the simplifications - but the first calculation needs 
specification of n as an integer. 

I believe that the “not-simplifying” is connected with the internal treatment as complex numbers. 


n n 


2 -n 


# 1 : (- 1 ) .(- 1 ) = (- 1 ) 

#2: n :e Integer 

n n 

#3: (-1) •(-!) = 1 


2 - n ! 

#4: = 2 

n! 

(n + 3)! 


#5: 


= (n + l)-(n + 2)-(n + 3) 


An interesting mathematical request: 


Hi; 

My objective is to discover the exact 3D points for creating Platonic solids. Since these are the only geometries 
that have equal length sides and equal angles between them, I thought it might be easiest to find a formula that 
enables me to calculate farthest points within a sphere where the lengths and angles are equal. Does such a for- 
mula exist? 

Jonathan Pike [jonathanpike55@yahoo.com] 

Hi; 

The Platonic solids nest perfectly one within another. This configuration, starting with the octahedron, is known 
as Metranon's Cube. I would like to build a mathematically perfect model of the same. Andrew Robbins was 
kind enough to point me to this page: 

http://mathworld.wolfram.com/PlatonicSolid.html 

which gives me a lot of pertinent data; however, the data for exactly nesting one solid in another seems lacking 
by half. It gives the formulas for determining the spherical radii of nested solids vis-a-vis their "dual" solids; 
however, those duals are matched icosahedron-dodecahedron, octahedron-cube and tetrahedron-tetrahedron (for 
a star tetrahedron). In the Metranon's Cube as described above, the tetrahedron is between the octahedron and the 
cube, so how does one figure out the "duals" or the radii of nested solids? This problem compounds itself with 
the matter of bridging the cube to the icosahedron. How do I resolve these matters? 

Jonathan 
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A TI-N spire problem (conflict?) with radian and degree 

I defined a straight throw in parameter form dependent on the angle. I set the degree mode in the gen- 
eral settings. Surprisingly I obtain a graph showing a negative slope (with an angle of 30°) in the ge- 
ometry window. The function table is correct - fitting to the degree settings. Working in rad mode the 
graph appears correctly. 

My conjecture is: the graph is based on the radiant although the equation gives correct function values 
which are given in the table. 

My question is: Is my conjecture right or am I doing something wrong? 

Many thanks and best regards 
Helmut 



The “trick” is: You have to set the Geometry Angle to Degree in the Application settings for 
the Graphs & Geometry Application,too. It is not sufficient doing this in the General Settings. 
Don’t ask me why one has to do this twice. Josef 
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Robert Setif s “autistic problem” 

Dear Josef, 

My function naff 3 is not correct. I can't master the command append (out, res 
sub [ ,...,] nor append (out, [res sub [ , . . . , ] . I would obtain 123 456 or 

[123], [456] or [ [123] , (456] ] for 123456; 

I would like to have 1 234 567 890 for 1234567890 or with the minimum of 
[] 12 345 678 901 for 12345678901. 

Best regards. 

Robert 

naff 3(1234567891, 3) = [[[1, 234, 567, 891]]] 

Von: Robert SETIF [robert.setif@gmail.com] 

Gesendet: Dienstag, 17. Mai 2011 12:41 
Betreff: naffndigs ! ! ! 

Dear Josef, 

Excuse me with my "autistic" problems. 

I would like a function that displays only the fractional part . (or the 
integer part and the rest below) . 

For 3.14... -> 14159 26535 89793 ... or 141 592 653 ... 

This is not urgent. 

Thank you very much. 

Best regards. 

Robert 

Dear Robert, 

in earlier times a used to program a lot using string operations. 

I attach my version of diffn.mth. 

The trick is to insert a blank (= “ “) between the groups of characters. 

Best regards 
Josef 


naffndigs(n, digs, res, out) := 

Prog 

ni= STRING(n) 
out := 11 " 

Loop 

#1: If DIM(n) < digs 

Prog 

out := APPEND(n, out) 

RETURN out 

out := APPENDC ", m[DIM(n) - digs + 1 DIM(n)] , out) 

n := n i [ 1 DIM(n) - digs] 

#2: naffndigs(1234567890 , 4) = 12 3456 7890 

#3: naffndigs(1234567890 , 3) = 1 234 567 890 

#4: naffndigs (1234567890, 4) = 12 3456 7890 

#5: naffndigs(1234567890, 6) = 1234 567890 

#6: naffndigs(1234567890, 5) = 12345 67890 

#7: naffndigs (1234567890 1234567890 , 3) = 12 345 678 901 234 567 890 

#8: naffndigs (1234567890 1234567890 , 2) = 12 34 56 78 90 12 34 56 78 90 

#9: naffndigs (1234567890 1234567890 , 1) =12345678901234567890 
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Rick Nungester wrote in the eDUG-Fomm: 

I am using Derive 5.06 for Windows, looking at exponentiation with negative base and/or power. "#i" 
is "i-hat" or SQRT(-l). Here are Derive results: 

-0.62 A 3.4 

Simplify: -29791*38440 A (l/5)/1250000 
Approximate: -0.1968483343 

Same base value, but complex representation: (-0.62 + 0*#i) A 3.4 
Simplify: 

-2979 1 *46 1 7605000000 A (l/10)/5000000+2979 1 *38440 A (l/5)/5000000 +/- 

#i*(559358543 18 121 35 14687266633 126576650487061 55 17773 13 1*SQR(5)/4768371 58203 1250000 
00000000000000000000000000000000000000000+5085077665283759224429696664787060458973 
2868343392 1/1 9073486328 125000000000000000000000000000000000000000000000) A ( 1/10) 

Approximate: -0.06082948062 +/- 0. 1 872 1 389 1 0*#i 

*** Why aren't the two answers the same? 

*** Are both answers correct? 

I thought complex exponentiation z A u is defined as e A (u*LN(z)), or as an algorithm, find LN(z) (a 
function with principal branch as part of its definition and undefined at z = 0), then find u*LN(z) (a 
function), then find e A (u*LN(z)) (a function). 

I also thought x A y is real if (x > 0) or (x = 0 and y > 0) or (x < 0 and y is integer), with 0 A 0 depending 
on who is doing the inteipretation. For example, 0 A 0 = 1 in Derive but = "undefined" on the HP-32SII 
calculator. 

*** Where else am I confused? 

The above is a typing error. It should be: (-0.62) A 3.4 
I got this reply via email: 

...here (-0.62 A 3.4) approximates to -0.06082948062 - 0.187213891(H(same version ofDERIVE). The 
absolute value of the above approximates to 0.1968483343, in accordance with your post. 

My guess is that some non-standard settings (about branches and cuts in the complex plane) are in 
effect on your system. Try starting up DERIVE with the default (factory) settings and repeat the calcu- 
lation. 

Dirk 


Thanks Dirk, that was it. I had Derive set to not prompt me if I wanted to start with factory 
settings, and have been accumulating non-factory changes for years. I don’t know how I got 
my posted results, but now don’t care. From now on I plan to always start with factory de- 
faults and only deviate from them intentionally, saving those commands as part of my saved 
.dfw files. 
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A Mathematical Model for Snail Shells (3) 

Piotr Trebisz, Germany 

In meinen ersten zwei Beitragen iiber die Trebisz-Spiralen vom Typ 1 und 2 bin ich ausfiihrlich auf die Kon- 
struktion von Schneckenhausern eingegangen die auf planaren bzw. konischen logarithmischen Spiralen basie- 
ren. Man konnte also meinen dass damit schon alles iiber beide Schneckenhaus-Typen gesagt wurde. 

Beiden Schneckenhaus-Typen ist gemein, dass wir das lokale mitdrehende Koordinatensystem mit Hilfe der 
frenetschen Formeln hergeleitet haben. In der Differentialgeometrie ist das die Standardmethode fur die Be- 
stimmung eines begleitenden Dreibeins einer Raumkurve. 

Ich werde aber in diesem und im nachsten Beitrag zeigen dass es zu jeder planaren und konischen Spirale ein 
korrespondierendes orthogonales krummliniges Koordinatensystem gibt, das sich aufganz natiirliche Weise aus 
der Spiralen-Gleichung ergibt und in dem das Selbstdurchdringungsproblem auf einfache Weise analytisch los- 
bar ist! 

In diesem Beitrag werde ich mich auf planare logarithmische Spiralen konzentrieren. Die daraus abgeleiteten 
Schneckenhauser nenne ich "Trebisz-Spiralen - Typ 3". 

In krummlinigen Koordinatensystemen sind die Koordinatenachsen keine Geraden sondern Raumkurven. Die 
wohl bekanntesten sind das spharische und das zylindrische Koordinatensystem. 

In my first contributions I treated snail shells based on planar and conical logarithmic spirals. Both 
types have in common that the local co-rotating system of coordinates was derived using the Frenet 
formulae (which is the standard method in differential geometry). 

In this contribution and in the next one as well I will demonstrate that there exists a corresponding 
curvilinear coordinate system to each planar and conical spiral, which can be easily derived from the 
spiral eguation. Then the interpenetrating problem can be solved analytically. 

Here I will concentrate on planar logarithmic spirals. The resulting snail shells will be called “Trebisz- 
Spirals - Type 3”. 

In a curvilinear coordinate system the coordinate axes are no straight lines but space curves. The best 
known of these systems are the spherical and the cylindrical coordinate system. 


Spharisches Koordinatensystem: 


Zylindrisches Koordinatensystem: 
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Worin liegt nun also die Motivation dafur, ein krummliniges Koordinatensystem zu verwenden? Dazu 
werfen wir am besten einen Blick auf den Langsschnitt durch ein Schneckenhaus vom Typ 1, das nach 
der Windungszahl parametrisiert ist. Die Spiralkurve, die als Basis fur das Schneckenhaus dient, ist 
hier in blau dargestellt. Die roten Linien sind Projektionen der sich beriihrenden Mantelkreise in der 
YZ-Ebene. 

What is the motivation for using a curvilinear coordinate system? Let’s have a look on the 
axial section of a type 1 snail shell, which is parameterized by the winding number. The spiral 
which serves as base for the shell is presented blue. The red segments are the projections of 
the osculating circles of the surface into the yz-plane. 

Beispiel: b = e und 0 < t < e 



Es zeigt sich deutlich, dass die Mantelkreise einen Kantenzug ergeben der wiederum eine logarithmi- 
sche Antispirale approximiert! Das erkennt man daran, dass der Winkel und das Langenverhaltnis 
zwischen zwei benachbarten Kanten konstant sind. Die beiden Spiralen schneiden sich in den Schnitt- 
punkten i mm er unter dem Winkel von 90°. Das ist auch nicht weiter iiberraschend, denn die Mantel- 
kreise stehen ja defmitionsgemaB immer senkrecht auf den Tangenten der Spiralkurve. 

Ein Schneckenhaus vom Typ 3 entspricht im Wesentlichen einem Schneckenhaus vom Typ 1. Als 
Geriist dient dieselbe planare logarithmische Spirale und der Radius des Mantelkreises wachst auch 
linear mit dem Radius bzw. mit der Lange der Spirale. Der Unterschied besteht nur darin, dass der 
Mantelkreis in einem krummlinigen lokalen Koordinatensystem defmiert ist und nicht in einem karte- 
sischen begleitenden Dreibein wie beim Schneckenhaus vom Typ 1. Aus dem Kantenzug im Langs- 
schnitt wird daher eine perfekte, stetige Antispirale. 

We can see that the segments approximate a logarithmic antispiral. This can be recognized 
by the fact that the angles and the ratios of consecutive segments is constant. Both spirals 
intersect with 90°. This is not really surprising because the surface circles are - by definition 
- always perpendicular to the tangents of the spiral. 

A type 3 snail shell basically corresponds with a type 1 shell. The same planar logarithmic 
spiral serves as scaffold and the radius of the surface increases iineariy with the radius and 
the length of the spiral respectively. The difference is that the surface circle is now defined in 
a curvilinear local coordinate system and not in a Cartesian three pod (type 1). The segment 
progression changes to a perfect continuous antispiral. 


D-N-L#84 


Piotr Trebisz: A Mathematical Model for Snail Shells (3) 


p 41 


Bei spi el : b = e und 0 < t < e 



Die Verwendung einer solchen Antispirale als lokale Y-Achse bietet aber nicht nur asthetische, son- 
dern auch ganz handfeste Vorteile. Die Beruhrungspunkte der Mantelkreise liegen hier alle auf einer 
stetigen Kurve, dadurch wird das Selbstdurchdringungsproblem analytisch losbar! 

Ausgangspunkt der Herleitung ist die eigentliche Spirale, die als Basis fur das Schneckenhaus dient, 
und ihre Tangentengleichung. Der Parameter t steht hierbei fur den Radius bzw. die Lange der Spirale, 
wahrend b die Basis des Logarithmus ist. 

Using such an „antispiral“ as local y-axis does not offer only aesthetic benefits but also con- 
crete ones. All the osculation points are lying on a continuous curve which makes the inter- 
penetrating problem solvable. 

Starting point of the derivation is the spiral which serves as base for the shell and its tangent 
equation. Parameter t describes the radius and the length of the spiral, respectively, and b is 
the base of the logarithm. 

#3: SPIRALE (b, t) := [0, t -COS (2 tt- L0G(t , b)), t ■ SIN (2 -tt- L0G(t , b))] 

d 

#4: SPI R AL E_TANG E NT E ( b , t) := — SPI RALE (b, t) 

dt 

Gesucht ist eine logarithmische Antispirale, die die Basisspirale immer senkrecht schneidet. 

We are looking for an antispiral which intersects the base spiral always orthogonally. 

#5: ANTI_SPIRALE(a, s) := [0, s -COS (2 -tt- L0G(s , a)), s ■ SIN (2 -tt- L0G(s , a))] 

d 

#6: ANT I _S P I RAL E_T ANG E NT E ( a , s) := — ANTI_SPIRALE(a, s) 

ds 

Es erweist sich hier als sehr hilfreich, dass planare logarithmische Spiralen bei t= 1 bzw. 5=1 immer 
durch den Punkt [0, 1, 0] gehen, unabhangig von ihrer Basis. 

It proves to be helpful that planar logarithmic spirals for t = 1 and s = 1 are passing the point 
[0, 1, 0] independent on their base. 

#7: SPIRALE(b , 1) = ANTI _SP I RAL E ( a , 1) = [0, 1, 0] 

Gesucht ist eine Basis a fur die Antispirale, so dass die beiden Tangenten in den Schnittpunkten im- 
mer senkrecht zu einander stehen. 

We try to find a base a for the antispiral such that both tangents in the intersection points are 
perpendicular. 
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BASIS(b) := (SOLUTIONS (SRI RALE_TANGENTE(b, 1) • ANTI_SPIRALE_TANGENTE(a, 1) = 0, a, Real)) 
# 8 : 1 

2 

#9: - 4-tt /LN(b) 

BASIS(b) = e 

Die Lange der Antispirale hangt linear vom Parameter 5 ab, da ihre Tangente an jedem Punkt eine 
konstante Lange aufweist. 

The length of the antispiral depends linearly on the parameter „s because its tangent is of 
constant length in every of its points. 

2 2 

./(LN(b) + 4 ■ 7T ) 

#10: | ANTI_SPIRALE_TANGENTE(BASIS(b) , s) | = 

2 -7T 

Da die Anti-Spirale also nach ihrer Lange parametrisiert ist, aber multipliziert um einen konstanten 
MaBstabsfaktor, ist sie als Y -Achse in einem krummlinigen Koordinatensystem gut geeignet. 

Als X-Achse benotigen wir einen Vektor, der senkrecht zu der Ebenen steht in der sich die Basisspira- 
le windet. Auf diese Weise stellen wir sicher, dass die lokale X-Achse immer se nk recht auf die Tan- 
genten der Antispirale und der Basisspirale steht. 

Die Gleichungen fur die Achsen des krummlinigen Koordinatensystems lauten daher: 

As the antispiral is parameterized wrt to its length but multiplied by a constant scaling factor, 
it is well suitable for serving as y-axis in a curvilinear coordinate system. 

As x-axis we need a vector perpendicular to the plane which contains the base spiral. So we 
can make sure that the local x-axis is always perpendicular to the tangents of base spiral and 
antispiral as well. 

The equations for the axes of the curvilinear coordinate system are: 

#11: X_ACHSE(b, x) := | ANTI_SPIRALE_TANGENTE(BASIS(b) , s) | ■ [x , 0, 0] 

#12: Y_ACHSE(b, y) := ANTI_SPIRALE(BASIS(b) , y) 


Um zu gewahrleisten dass der MaBstab in Richtung der X-Achse genauso groB ist wie in Richtung der 
Y -Achse, wird die X-Achse mit dem konstanten MaBstabsfaktor der Y -Achse multipliziert. 

Zudem ist dieses krummlinige Koordinatensystem orthogonal, das heiBt, dass die partiellen Ablei- 
tungsvektoren immer senkrecht zu einander stehen. 

To make sure that the scaling in direction of the x-axis is equal to the scaling in direction of 
the y-axis the x-axis is multiplied by the constant scaling factor of the y-axis. 


Additionally this curvilinear coordinate system is orthogonal, i.e. the vectors of the partial 
derivatives are perpendicular. 


#18: 


fd ) 

— OLACHSEfb, x) + Y_ACHSE(b, y)) 
t dx v 


d 

• — (X_ACHSE(b, x) + Y_ACHSE(b, y)) = 0 
dy 


Der Mantelkreis wird im krummlinigen Koordinatensystem defmiert. Sein Radius soli linear mit dem 
Radius bzw. mit der Lange der Basis-Spirale wachsen, m ist dabei der Linearfaktor. 


The surface circle is defined in the curvilinear coordinate system. Its radius shall grow linear 
with the radius and the length of the base spiral, respectively, with m as linear factor. 






D-N-L#84 


Piotr Trebisz: A Mathematical Model for Snail Shells (3) 


p 43 


#14: KREIS(b, m, p) := X_ACHSE(b, m- COS (<(>)) + Y_ACHSE(b, 1 + m-SINOtO) 

Damit haben wir die Gleichung des Mantelkreises fur t= 1 dessen Zentrum im Punkt [0, 1, 0] liegt. 
Das lasst sich zeigen indem man den Radius auf 0 reduziert. 

The result is the equation of the surface circle fort = 1 with its centre in point [0, 1, 0], This 
can be shown by reducing the radius to 0. 

#15: KREIS(b, 0, <(j) = [0, 1, 0] 

Um daraus die Mantelflache des Schneckenhauses zu konstruieren, muss man den Mantelkreis um den 
Winkel der Basisspirale rotieren und mit ihrem Radius skalieren. 

For constructing the surface of the snail shell one has to rotate the surface circle by the angle 
of the base spiral and scale it by its radius. 

#16: MANTEL (b , m, t, <f0 := t ■ R0TATE_X(2 -jt- L 0G(t , b)) • KREIS(b, m, <f0 
#17: MANTEL (b , 0, t, <f0 = SPIRALE(b, t) 

Die Beruhrungspunkte der Mantelkreise liegen nach jeder Windung auf der Antispiralen. Dieser Um- 
stand ermoglicht es uns, das Selbstdurchdringungsproblem analytisch zu losen! 

Wir wissen, dass sich die Zentren der Mantelkreise nach jeder Windung in den Schnittpunkten der 
Basisspirale mit der Antispirale befmden. 

Es miissen also zuerst die Schnittpunkte der Spiralen berechnet werden. 

Sowohl die Basisspirale als auch die Antispirale sind nach ihrem Radius parametrisiert. 

The osculation points of the surface circles are located on the antispirals after every turn. 

This fact makes possible solving the interpenetration problem analytically! 

We know that the centres of the surface circles after every turn are located in the intersection 
points of the base spiral and the antispiral. 

So we have to calculate the intersection points of the spirals first. 

The base spiral and the antispiral are both parameterized wrt to their radius. 

#18: | SPIRALE(b , t) | = t a | ANTI_SPIRALE(BASIS(b) , s) | = s 

Wir konnen deshalb fur beide Spiralen denselben Parameter t verwenden, was uns die Berechnung der 
Schnittpunkte sehr vereinfacht. 

So we can use the same parameter t for both spirals which makes calculation of the intersec- 
tion points much easier. 

SOLUTIONS((SPIRALE(b, t)) = (ANTI_SPIRALE(BASIS(b) , t)) a 

#19 : 2 2 

SOLVE((SPIRALE(b , t)) = (ANTI_SPIRALE(BASIS(b) , t)) , t, Real), t, Real) 

3 3 


2 2 2 2 2 2 
#20: - 4 ■ 7T /(LN(b) + 4-7T ) A-n /(LN(b) + 4-7T ) 

_0, 1, b , b 

Die Basisspirale und die Antispirale haben aber nicht nur vier, sondem unendlich viele Schnittpunkte. 
Anhand der vier speziellen Losungen, die uns DERIVE liefert, konnen wir eine allgemeine Losung fur 
die Schnittpunkte angeben. k ist hier eine beliebige ganze Zahl. 
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The base spiral and the antispiral don’t have only four intersection points but infinitely many. 
Inspecting the four special solutions delivered by DERIVE we can indicate a general solution 
for the intersection points, k is an arbitrary integer. 

2 2 2 
#21: 4 ■ 7T ■ k/ CLNCb) + 4-rr ) 

SCHNITT (b , k) := b 

Damit konnen wir den Linearfaktor m so bestimmen, dass sich die Mantelflache nach jeder Windung 
exakt beriihrt ohne sich selbst zu durchdringen. 

Hence we can define the linear factor m in such a way that the surface is osculating itself 
after each turn without interpenetrating. 


M_YZ(b) := (SOLUTIONS(SCHNITT(b, k) ■ (1 + m) = SCHNITT (b , k + 1) ■ (1 - m) , m, Real)) 
# 22 : 1 

2 2 2 

4-rr /CLNCb) + 4-rr ) 
b - 1 

#23: M_YZCb) = 

2 2 2 

4-rr /CLNCb) + 4-rr ) 
b +1 

Aus der Differenz zweier aufeinander folgender Schnittpunkte konnen wir auBerdem direkt ablesen, 
wie groB eine Windung beim Schneckenhaus vom Typ 3 ist, also die Umdrehungszahl nach der sich 
die Mantelflache selbst beriihrt. 

The difference of two consecutive intersection points enables reading off how large one turn 
of a type 3 snail shell is, i.e. the winding number after which the surface is osculating itself. 

#24: WINDUNG Cb) := LOGCSCHNITT Cb , k + 1) , b) - LOGCSCHNITT Cb , k) , b) 

2 

4-77 

#25: WINDUNG Cb) = 

2 2 
LNCb) + 4 ■ 7T 

Und damit sind wir fertig! Man muss nur das passende m in die Mantelflache einsetzen und erhalt so 
ein Schneckenhaus vom Typ 3, das nach dem Radius bzw. nach der Lange parametrisiert ist. 

Now it is done, we are ready! We have to substitute the right m into the surface expression in 
order to obtain a type 3 snail shell which is parameterized wrt to its length. 

#26: SchneckenHausYZ_RCb , t, tp) := MANTELCb, M_YZCb), t, (p) 

Will man das Schneckenhaus nicht nach dem Radius, sondern nach der Zahl der Umdrehungen 6 pa- 
rametrisieren, dann ersetzt man t durch b e . Der Radius wachst dann nicht mehr linear zu t sondern 
zu b°. 

If one does not want parameterize the snail shell wrt the radius but wrt to the number of turns 
0 then one has to replace t by b 0 . Then the radius will grow linear to b 0 . 

e e 

#2S: SPIRALECb, b ) = b 

9 

#29: SchneckenHausYZJJCb, 6, (p) := MANTELCb, M_YZCb), b , tp) 
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Da wir wissen welcher Umdrehungszahl eine Windung entspricht, konnen wir auch ein Schnecken- 
haus definieren das nach der Anzahl der Windungen parametrisiert ist. 


As we know which rotation number corresponds with one turn we can define a snail shell 
parameterized by the number of turns. 


#31: SchneckenHausYZ_W(b , 8, tp) := SchneckenHausYZJJ(b , 8- WINDUNG(b) , tp) 

Und damit die ganze Rechnerei nicht umsonst war werden wir belohnt mit einigen schonen Schne- 
ckenhausem. 

Finally we are rewarded for all our calculations with some pretty snail shells. 


Basis 2: 


t 

#40: MANTEL (2, M_YZ(2), 2 , s) 


(-71 < S 71, -4 < t < 1 ) 



Base 6: 


t 

#41: MANTEL (6, M_YZ(6), 6 , s) 


(-71 < S 71, -2 < t < 1 ) 
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Base e with exact one turn: 


t 

#42: MANTEL (e, M_YZ(e) , e , s) 


(-7t < S 71, 0 < t < 1 ) 



You can study the influence of b by introducing a slider bar for b, Josef. 



This is another spiral. 

Do you recognize which one? Josef 
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Erik van Lantschoot: TRJGO for TI-NspireCAS 
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Erik van Lantschoot: TRIGO for TI-NspireCAS 
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There is an auxiliary function wstw hidden in the program: 

Define wstw(g)= 

Func 

:© wstw=where is what? g=list of boolean expression, e.g. where do I find 
:© a side and no angle? then: g:={not ks and kw}. (ks = keine Seite = no side) 
:© The indices of the positions of "true"in g are collected in list 1st. 

:Local h,lst: lst:={} 

:Forh,1,3: If g[h]=true Then: lst:=augment(lst,{h}): Endlf: EndFor 

: Return 1st 

:EndFunc 


I wrote above that it should be no problem to transfer TRIGO to the other Tl- 
calculators. 

I tried for the Voyage 200. The most important difference is that assigning values 
needs a -> on the Voyage (TI-92) instead of := with Nspire. 

See some screen shots from the V200 treating the same problems as before with 
Tl-Nspire CAS. 

Josef 
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A Lucky Charm for the DU6 
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Roland Schroder's gift for the DUS Members: 

Roland sent a DERIVE made Four-Leaf Clover as a Lucky Charm 2012 for all DUS 
Members. 




f 

n 

> 



HOD (5 , 103). 2- rr 


VECTOR 


1, 

, n, 0, 103 


■v 

103 

J 


??? 


I wrote back: 

I did an extra work looking for the envelope of the family of lines. It fits exactly. 

Thanks for the challenge! 

(I don't present the formula here and leave it for our readers. It will be presented in DNL#85, Josef) 

Roland answered (in Ger man and I translate) 

Dear Josef, 

I find it really striking that your set of points describes the four-leave clove. As I had 
never to do geometry I don't have any clue why this is the case and how to find such a 
formula. 

I found my formula during a longer “expedition' 1 which started with the question: What 
will happen when representing sequences of numbers on a circle instead on a number 
ray? As you might guess this results in a modular reduction of the sequences. And this 
leads us to in number theory. 

You can find the sequence of powers of two and three presented on a circle in one of my 
papers of DERIVE applications "Light in the coffee cup" . Some of these papers did you 
publish in earlier DNLs. 

The reason for that lies in the depth of number theory. There is the concept of the 
“primitive root": the constant base of my sequence of powers must be primitive root of 
the prime number which indicates the number of points on the circle (the module). You 
can find a website which gives the respective primitive roots. You will notice that 
square numbers cannot be primitive roots (which can easily be proved). So 4 cannot be 
primitive root of the number of points on the circle. It follows that a clover leaf with 4 
-1 = 3 leaves is not easy to create. 

Best regards 
Roland 
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More Fractals produced with Fractalizer 
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Download the freeware fractalizer from www.fractalizer.de , Josef 







